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We present numerical results for the zero temperature persistent currents carried by interacting 
spinless electrons in disordered one dimensional continuous rings. The disorder potential is described 
by a collection of (5-functions at random locations and strengths. The calculations are performed 
by a self-consistent Hartree-Fock (H-F) approximation. Because the H-F approximation retains the 
concept of single-electron levels, we compare the statistics of energy levels of noninteracting electrons 
with those of interacting electrons as well as of the level persistent currents. We find that the e-e 
interactions alters the levels and samples persistent currents and introduces a preffered diamagnetic 
current direction. In contrast to the analogous calculations that recently appeared in the literature 
for interacting spinless electrons in the presence of moderate disorder in tight-binding models we 
find no suppression of the persistent currents due to the e-e interactions. 



I. INTRODUCTION 



The prediction of Biittiker, Imry and Landauerl that a 
persistent current (PC)Bfl can be observed in a disordered 
normal metal mesoscopic ring threaded by a magnetic 
flux, although this ring has a finite resistance due to ( 
tic scattering, has attracted much recent attentior 
The elastic scattering, due to disorder, reduces the am- 
plitude of this current. The detection of the effect in 
three different experiments, by measuring the magnetic 
signal of the PC, has stimulated a great deal of theoret- 
ical interest in particular because of the large magnetic 
response that was measured. The measured amplitude of 
the PC, although some of the samples are clearly in the 
disordered diffusive regime, is of the same magnitude as 
that calculated for electrons in a clean ring. 

The simultaneous measurement^ of the magnetic re- 
sponse of 10^ metallic Cu rings yields PC amplitude 
which is one order of magnitude larger than predicted 
by single-electron calculations which take into account 
the elastic scattering due to the disorder. This rised the 
question of the source of such large magnetic response 
of an ensemble of rings which are in the diffusive regime 
of disorder. The experimental signals of three single Au 
ringsllil of different geometries were found to be up to two 
orders of magnitudes larger than the prediction of such 
single electron theory. On the other hand, the responseEj 
of a 2D semiconducting ring in the ballistic regime of dis- 
order and at a very low electronic density agrees well with 
the theoretical predictions of the single electron theory. 
These experimental findings rise the question whether 
the discrepancies between single electron theory and the 
observed behavior of the PC can be understood as a con- 
sequence of neglecting the e-e interactions in the pres- 
ence of disorder. This question is part of a more general 
problem of the influence of electronic interactions in the 
presence of disorder on the properties of mesoscopic sys- 
tems which is in the center of many recent studies, such 



as polarizabilitji fluctuationsEl, quasiparticle life times in 
quantum dotsllj and the metal-insulator transitiorJi3. 

There have been many attempts to investigate the influ- 
ence of e-e interactions on the PC. Analytical calcula^ 
tions were mainly restricted to perturbation theoriest^ 
and tO|-^enormalization treatments of Luttinger liquid 
modelsim, which have shown that some enhancement of 
the PC due to the e-e interaction is possible. Exact 
numerical diagonaJization studies of spipJess interacting 
electrons for lEoO as well as for 2DE3 tight-binding 
(TB) models were performed. In ID systems for long 
range Coulomb interaction the PC is weakly enhanced 
for strong disorder (i.e., for the localized regime) and 
medium strength of interactions. But for weak disorder 
the current is suppressed at any strength of the interac- 
tion and filling factortS. Also for any strength of short 
range e-e interactions^ven far from half filling, the PC 
has shown a decreaseEj. On the other hand, for 2D sys- 
tems tJbe PC is significantly enhanced in the diffusive 
regimec2l. Although these studies are useful in point- 
ing towards the general influence of the e-e interactions 
they are restricted, due to exponential expansion of the 
Hilbert space with sample size, to samples of small num- 
ber of sites (about six to twenty) . This complex behavior 
has its origin in the interplay between two different phys- 
ical effects. While the disorder creates fluctuations in the 
electronic density leading eventually to localization, the 
e-e interaction is expected to reinforce the tendency for 
local charge neutrality which counteracts the influence 
of disorder. In the H-F approximation this is realized 
by the effective potential which is expected to lead to 
an increase in the PC amplitude. On the other hand 
the Mott-Hubbard transitiorO reduces the current am- 
plitude because the e-e interaction opens a gap at the 
middle of the conduction band leading to an insulating 
-not metallic- behavior of a half flUed band. 

It is interesting to examine the PC behavior in the 
H-F approximation because exact analytical considera- 
tion of interaction in the presence of disorder is a very 
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difficult task. Even perturbative calculations are not 
conclusiveE^. On the other hand exact numerical diag- 
onalization is limited to small number of sites and elec- 
trons. But using the H-F approximation has the advan- 
tage of the possibility to handle much larger systems. It 
has also the possibility to differentiate between different 
influences of the e-e interactions which may lead to an 
enhancement of the PC. Further it retains the concept 
of single electron levels in a meaningful way. This en- 
ables one to employ the considerable existing knowledge 
on statistical properties of single electron energy levels in 
disordered systemsE3 for studying the influence of e-e in- 
teractions. One is also able to study the behavior of the 
current for different regimes of disorder since as the sin- 
gle electron energy becomes higher the effective disorder 
for these energies is weaker. 

Kato and Yoshiokall were the first to deal with large ID 
TB system of 100 sites and about half filling in the frame 
work of the self consistent H-F approximation for long 
range interaction. They have found that for any value of 
disorder the electronic interactions further suppress the 
PC and do not counteract the effect of disorder. This 
difference with the results of ID exact diagonalization is 
due to the fact that the HF approximation for spinless 
electrons does not smear the density fiuctuations which 
is essential to obtain an enhancement of the curi:ent in 
the localized regime. For spinless electrons in twco and 
threeE3 dimensional TB rings the HF calculations show 
some enhancement of the average current for Coulomb 
interactions in the diffusive regime. But if the spin is 
taken into account the HF approximation shows a large 
enhancement of the PCcB. The reason is that including 
the spin degrees of freedom adds a dominant contribution 
to the direct term. In fact, this pffptjt can be seen even 
for first order perturbation theorycZlH. Thus, for the case 
where spin is taken into account, the HF approximation 
does describe the appearance of local charge neutrality 
due to interactions. 

The main goal of this paper is to compare the PC behav- 
ior of interacting spinless electrons in continuous disor- 
dered rings with that behsudoi-iixi TB rings that recently 
appeared in the literatureHll^O. For these TB rings in 
the moderate regime of disorder it was found that e-e in- 
teraction further suppresses the PC because of the Mott- 
Hubbard transition which does not seem to be relevant 
for electrons in continuous rings. The difference in the 
behavior of the PC is a measure for the importance of the 
difference between continuous and discrete models. The 
difference between a model with no discrete symmetries 
and the TB model has been suggested in Ref. and 
here we shall examine this difference for a well defined 
Hamiltonian. 

The rest of the paper is organized as follows: in sec- 
tion II we present the model for independent electrons 
and interacting electrons within the self consistent HF 
approximation. In section III we discuss the numerical 



results, and in section IV we derive our conclusions. 



II. MODEL OF ELECTRONS IN DISORDERED 
ID CONTINUOUS RINGS. 

A. Non interacting electrons 

Let R be the radius of the continuous ID ring, the 
coordinate along the ring, me is the electron mass, e its 
charge, (j) the flux threading the ring and (j)o = ^ the 
flux quantum. For completeness and further reference 
let us restate some basic facts. Defining the disorder 
potential as Vdisid) the Schrodinger equation (in energy 
units = 1 ) for non- interacting electrons is: 

where the boundary conditions (BC) are periodic 

^-(0 + 27r) = *(6l). (2) 

The substitutionEI ^{0) fi V'(^')e~*^^ yields a wave 
equation of the Bloch typelil 

- ^^^(^) + yd^si0)m = mo), (3) 

which now has the modified BC 

ijie + 27r) = exp (i27r^)?/'(6'). (4) 

An exact eigenstate of energy Ei carries PC ii 

(where numerically we use dimensionless energy and di- 
mensionless current —-^M-)- We stress that \E'i(0) are 

the electronic wave functions which are exactly periodic 
for any value of (j). This fact will be important when 
we consider interactions. For Vdis{0) = the electronic 
fmictions arc ^i{9) = exp {il9) with / = 0, ±1, ±2, ... 

and the spectrum is Ei — ^xf with xi = I + For our 
model we have chosen 

Vd^s{0) = J2xAs-0,) (6) 

where the locations < 9j < 2tt and strength — A < < 
A of the individual j-th scatterer are random with the ap- 
propriate probabilities: Pg{9j) = ^ and P^{Xj) = 
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Ns is the total number of scutterers in the ring. Vdis{d) 
produces characteristics of-disordered samples as was dis- 
cussed by Imry and ShirencJ while truncating the Hilbert 
space. For non-interacting electrons one can exactly find 
the spectrum and eigen functions by means of a numer- 
ical transfer matrix technique. A transfer matrix that 
"propagates" the solution of eq.(||) from the left of the 
j-th scatterer to the left of the next j + 1-st scatterer is 
given by: 



(1 + ^)e''- 



-ix{9j 



(1 



(7) 



where the exponents contain the effect of free propaga- 
tion of waves between the scutterers. If we denote as 
a two component vector in which the upper component 
corresponds to the amplitude of the forward propagat- 
ing wave e"^ and the lower component to e^"^ then 

The concepts of localization length ^ and the Landauer 
conductance, which are very useful in the study of open 
systems, can be defined by comparing eq.(j^) with the 
general form of a transfer matrix: 



1/t* 
-ir/t) 



-{r/ty 

l/t 



(8) 



{t and r the complex transmission and reflection coeffi- 
cients obeying |rp -I- |ip = 1). One can now analytically 
find ^ as apfunction of energy and also find the appropri- 
ate scaling^ behavior of the fluctuations of the Landauer 
I 1^ 

resistance p = ^ (or conductance g = 1/p), where Ng 
plays the role of the length scale of the conductoril. For 

2 

weak disorder, where pj = fj « 1, and small Ns one 

expects an Ohm law behavior, i.e., disorder averaging 
should give < p >^< pi > Ns (< pi > is the aver- 
age Landauer resistance of an individual scatterer). One 
finds 



<p>^ 



A2 



(9) 



as long as < p > is of order unity or smaller. A com- 
parison of numerical results with the above equation is 
presented in Fig.|l|a. 

Because of localization, for strong disorder or at any large 
enough Ns (where 1 << < p >) we expect < In(l-l-p) >~ 
jNs ■ In our case we analytically calculate ^ to be 



=< ln(l + pi) >= ln(l + A^) - 2 + ^ tan"! A 



(10) 



where A = ^ and pi is the resistance of an individual 
scatterer. This result is confirmed by numerical calcula- 
tions presented in Figj^ where the slope is the inverse 
localization length. 



Let us now return to the discussion of PC in the closed 
ring. The discrete spectrum and eigen functions of the 
ring can be found using the definition of the transfer ma- 
trix (eq. (0)) and the BC (eq. (^)) as an eigenvalue 
equation: 



e *OVi. 



(11) 



It is enough to deal with the real (or imaginary) part 
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FIG. 1. (a) Scaling behavior of the disorder average Lan- 
dauer resistance < p > as a function of the number of scatter- 
ers Ns for an open chain of scatterers at weak disorder (A = 5, 
E = 5000, see text), (b) Scaling behavior of < ln(l + p) > 
of an open chain at strong disorder (A = 20, E — 200). The 
slope defines the inverse localization length i. The empty 
symbols in (a) and (b) are our analytical estimates by eqs.(0) 
and (^) respectively, the full symbols are results of numerical 
simulations of 200 realizations. 
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of the last eigenvalue conditions in order to obtain the 
eigen values because the other part is automatically ful- 
filled. The case of A^^ = 1 at any strength of A can easily 
be solved. For 1 << TV^ and appreciable strength of Xj 
the current can be calculated numerically and is found 
to be suppressed by many orders of magnitudes due to 
localization. 

For a given disorder the total current0 in an isolated 
(canonical) sample is the sum of the currents of all the 
occupied levels in the ring, which for electrons at zero 
temperature corresponds to: 



We 



1=1 



(12) 



For a sample connected to a reservoir (grand canonical) 
the sum is over all levels having energies lower than a 
given chemical potential fi: 



direct and exchange terms will cancel each other. Nev- 
ertheless, for numerical convenience, we treat each term 
separately and check that this cutoff has no influence on 
the results. 

In order to approximate the integro-differential equa- 
tion (Eq.(14)) by an ordinary differential equation we 



would like to approximate the exchange term by a term 
which can be represented as an effective potential. This 
may be achieved by using the almost closure relation, 

E/='i ^rC^O^K^*) ~ ^i^' in the Fock term, which 

is reasonable when Nf, ^ 1. The last sum has a finite 
width which is of order This means that the main con- 
tribution to the Fock term integral comes from 0' ^ 6. 
Therefore one can exclude the unknown function tpiO') 
from the integpand and replace it by ipid) as a multi- 
plicative factoiEa. Thus in this approximation the Fock 



term is replaced by ■ 



ro Jo ^(e-0,)2+^2 



de'il^{e) and 



eq.(|l4|) becomes an ordinary Schrodinger equation: 



(13) 



1 ^2 



(15) 



B. Interacting electrons 

We include e — e interaction within the self consistent HF 
approximation. In analogy to the notations in eqs.(|l|) 
- (^ our equation for the new HF single particle wave 
function 'ip{9) reads: 

I R r -We f2^ de' jjQ\ 

with the BC of eq.(^). vq = eh^/niee'^, and e is the 
dielectric constant, to be distinguished from the cutoff 
e. The first term in the above square brackets takes 
into account the effect of the neutralizing background 
charge. This results in a constant contribution to the 
potential which will not affect the PC. The second term 
corresponds to the Hartree (direct) term and the third 
contribution represents the Fock (exchange) term. The 
square distance between the particles is approximated 
by (6 - e'f = min[|6' - 6»'|2, (27r -\e- e'\f]. A cutoff 
was introduced in order to make the contribution of each 
term finite. If one ignores the background term there is 
no need to introduce a cutoff since the divergence of the 



Here V"e//(6')V'(£) contains the three terms in the square 
brackets of eq.(|3) with the approximated Fock term. 
The resuh of Y.i=i '^*i{d')'^i{0) for a typical realization 
is presented in Fig.|[ It is clearly seen that, due to the 
sharpness of the almost closure relation, approximating 
eq.(p^) by eq.(|l5|) is reasonable. The approximation is 
reasonable for all electrons. At low energies the wave 
length of the considered electron is much larger than the 
width ^ Xp (where Af is the Fermi wave length) and 
therefore the approximation is reasonable at low energies. 
For an electron near the Fermi energy E » Veff{0) 
(where Veff{9) is defined to vanish at its minimum and 
is of the order of A) the error as a consequence of the 
approximation can not create any significant change in 
the current of such electron. 

We solve eq.([l5|) self-consistently. For the n-th iteration 
the single electron wave functions from the previous iter- 
ation are taken and used to construct the effective poten- 
tial Veff{0) for the current iteration. For the first itera- 
tion (n = 1) we define Veff{0) by the non-interacting so- 
lutions, i.e., equivalent to setting Veff{9) = at iteration 
n — 0. The smooth Veff{0) is approximated by a con- 
stant potential in the intervals between the external im- 
purities. If Ns is small a finer subdivision of the intervals 
may be necessary. For the k-th interval 9k < < dk+i the 
effective potential is defined as Veff{9) = Veff{9k)- The 
orthogonal solutions, which vary from one interval to the 
other, are given by e"^^^^ where Zk = ^2 (14// (6*) — E) 
is either real or pure imaginary. As one is free to add a 
global constant to the energy we define V^j / (9) to be zero 
at its minimum. For that interval the general solution 
'4'{9) is a linear combination of exp[±ix9) with some real 
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FIG. 2. Presentation of the almost closure relation: 
E^Il ~ S{e' - e), see eqs.(0) and (|l|), for one 

ring of the disordered ensemble at convergence of the self con- 
sistency {Ns = 30 random scatterers, A = 14, A'^e = 42 elec- 
trons). The real part of the sum is concentrated at 9 6' . 
The imaginary part of the sum is vanishingly small. This fig- 
ure justifies the approximation of the HF integro differential 
eq.(|l^ by eq.(^^ that we self consistently solved. Note that 
the electronic wave functions are 27r periodic for any 

flux and so is the effective potential. 



positive X that defines the energy E-.= jl. Propagating 
■0(6') along the ring circumferenceE^, with an appropriate 



matching at the boundaries of all intervals, asd applying 
the BC of eq.(^ yield the discrete spectrumEll and eigen 
functions of eq.(nq) which can now be normalized. 



We shall now discuss the choice of parameters for the nu- 
merical solution of Eq. (^5|) . The number of electrons in 
a given realization of disorder was set by requiring that 
the energy of the highest occupied state I will be smaller 
than 312.5, i.e., xi < 25, for any value of flux. This re- 
sults in an average filling of 42 electrons per realization. 
This requirement was chosen in order to balance between 
computation time limitations and maintaining enough 
electrons for the validity of the almost closure relation. 
This choice allows the study of different energy regimes 
such as the almost ballistic (higher energies), moderate 
(intermediate energies) and the strong (lower energies) 
disorder limit. The strength of the interaction term is 
determined by the ratio — . For metals this ratio for a 

J ro 

reasonable experimental setting (say R = 0.1/zm) will be 
about 10^, while for semiconductors (which have small 
effective masses and large dielectric constants) the ratio 
is of order of one. Since we can treat only a limited num- 
ber of electrons, we are really considering the low den- 
sity limit of the problem, which corresponds better to 
the semiconductor case and therefore we chose — = 1. 
A comparison between the TB model parameters and 
our estimation of the strength of interaction in semicon- 
ducting rings can be achieved by recognizing that for the 
typical densities of semiconducting devices the ratio of 
the electrostatic interaction energy U to the kineic en- 
ergy t (the hopping matrix element) is of order of unity^ 
for semiconductor. The cutoff was chosen as e = 10~^. 
It was checked that enlarging the cutoff up to 10"'^ does 
not change the results within the limits of accuracy. As 
a convergence condition for the n-th iteration at a given 



flux we required: 'Zi'Ji I ^T' - l< 0.00025. For a 

typical realization this condition was fulfilled within ten 
iterations. 

In order to calculate the PC one can use two different 
approaches. The first is the derivative of the many par- 
ticle ground state energy as function of the fiux and the 
second is a direct application of the current operator on 
the many particle wave function. As we have confirmed 
both methods should give the same current, i.e.. 



„(") 



„("-!) 



d 



< «'^|H|*„ >=< '^Ml'^v > 



(16) 



Here, ly and are the H-F ground PC and H-F ground 
state, H is the exact Hamiltonian. In the calculation of 
the curijeat it is sometimes convenient to use the following 
relationEa 



(17) 



where Ei are the eigenvalues of eq.(|l5|). It is worth- 
while to note that in the H-F approximation the PC 
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is not simply the derivative of the sum of the H- 
F single particle eigenvalues Ei. The single level 
current for a particular level is defined as ii = 



< 



> 



We shall note that 



the PC and the eigen states of a clean ring are completely 
unaffected by the e-e interaction. 



III. NUMERICAL RESULTS 



In this section we shall describe the results of our self 
consistent numerical solution of Eq. (|l5|). We have con- 
sidered realizations with Ng — 30 randomly placed delta 
scatterers of strength A = 14 (see eq.(|^)) for all the nu- 
merical calculations. The number of electrons occupying 
a specific realization was defined by the noninteracting 
problem as was explained in the previous sub-section. For 
a typical realization was of order of 42±6. T4/ / (9) was 
calculated as a staircase potential at 2Ns points (see the 
previous section for details) and was found to be a rea- 
sonable approximation of the smooth effective potential. 
Our ensemble of disordered rings contained 150 different 
realizations. In the following we present the results of our 
study of the effect of the c-e interactions on (A) the level 
spacing, (B) the PC distributions, and (C) the canonical 
average and the single ring PC. 



A. Spectrum 

As mentioned in the introduction, one of the main advan- 
tages of the H-F approximation is that although one is 
dealing with a many particle problem the single electron 
levels still characterize the system. Thus, the immediate 
question arises whether e-e interactions imposes a tran- 
sition on the level statistics. Because in the H-F approx- 
imation the system is still characterized by a single elec- 
tron spectrum of an effective Hamiltonian (see eq.([l5|)), 
we were trying to answer this question by comparing the 
level spacing statistics of the eigenvalues of eq.(||) with 
that of eq.(^5|) when the condition for self consistency 
is satisfied. The level spacing statistics are considered 
for different regimes of disorder: (i) the extended regime 
(levels within a high energy strip) which corresponds to 
Wigner statistics for the noninteracting electrons; (ii) the 
localized regime (levels within a low energy strip) corre- 
sponding to the Poisson statistics for the noninteracting 
electrons. For these low energies the levels are clearly lo- 
calized because the calculated levels currents were found 
to be exponentially small and by drawing \tjji{9)\'^ we sim- 
ply visualized the localization. Fig. ^ shows the results for 
both regimes. There is no transition in the level spacing 
due to the e-e interactions. The same analysis was carried 
out for an intermediate regime (levels in a regime between 
(i) and (ii)) also did not show any difference between the 
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FIG. 3. Comparison of level spacing statistics P(s) 
of noninteracting (empty symbols) and interacting (bold 
symbols) electrons within the H-F approximation, (a) At 
high energy strip P{s) is Wigner-Dyson N{l3)s'^ exp{ai3s'^). 
/3 = 1, ai = 7r/4, N{1) = 7r/2 for zero flux (dashed), 
l3 = 2, 02 = 4/7r N{2) = 32/tv) for half quantum of flux 
(long-dashed), (b) At low energy strip the states are local- 
ized also when including interaction as can be seen from the 
statistics which remains Poissonian. Because the change from 
(a) to (b) is continuous as a function of energy and P{s) is 
unchanged when including interaction there is no sign of in- 
teraction induced insulator to metal transition. 



statistics of the interacting and noninteracting electrons, 
although in both of these cases the expected transition 
from Poisson to Wigner statistics as a function of en- 
ergy was found. Thus, the level statistics do not show 
any sign of an insulator to metal transition, which may 
be expected if the interactions induce a straight forward 
change in the localization properties of the H-F levels. 
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Therefore, if there would be any enhancement of the PC 
(of this spinless model) it can not be explained by an 
interaction induced insulator to metal transition. 



B. Distributions of PC 



In order to understand the effect of e-e interactions on 
the canonical PC we plot the distribution of the canon- 
ical persistent currents for an ensemble of realizations 
with non-interacting electron as well as for the same re- 
alizations when interactions are taken into account. As 
can be clearly seen in Fig.^ the distribution of sample 
current is shifted towards more negative values, i.e. the 
PC becomes diamagnetic. Nevertheless, as a result of the 
interactions, individual realizations can change the cur- 
rent size and current direction in both the paramagnetic 
and diamagnetic directions. Unfortunately, because our 
ensemble is small it is meaningless to define the enhance- 
ment factor of the typical ^ {P) PC since it is strongly 
controlled by rare events in the tails. This shift of the 
distribution is similar to the situation seen in an exact 
diagorwjjzation study of small 2D systems iurthe diffusive 
regimec2l, although the direction is oppositeEZl. 

To get an information on the single electronic levels 
we have also compared the probability distributions of 
the single level PC for noninteracting and interacting 
electrons. One should separately analyze the different 
regimes of disorder previously defined. In Fig.||a the 
distributions of single level PC are given for two dif- 
ferent high energy strips in the extended regime. It is 
evident that due to the interaction the single level PCs 
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FIG. 4. Probabilities of sample PC for interacting (bold 
symbols) and noninteracting (empty symbols) electrons in 
the same canonical ensemble of 150 disordered rings at 
4- = 0.225. 
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are slightly shifted toward more negative values of the 
currents. The dip at the center of the probability distri- 
bution of the levels belonging to the highest energy strip 
is due to the fact that the electrons are almost ballistic 
and have a greater probability to carry large positive or 
negative current. In the lower energy strip the distribu- 
tion follows ci-Gaussian shape as was found by Simons 
and AltshulerO for levels in the diffusive regime. In con- 
trast to the high energy regime the effect of the interac- 
tions on electrons in the localized regime is to suppress 
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FIG. 5. Probabilities for level PC of interacting (bold 
symbols) and noninteracting (empty symbols) electrons in 
the same canonical ensemble of 150 disordered rings at 
^ = 0.225. (a) For extended levels within two high energy 
strips. At the highest strip the states are close to ballistic 
and have greater probability to carry large positive or nega- 
tive current. At the lower strip the probability is compared 
with a Gaussian (average and variance from numerical data) 
as predicted by Simons and Altshuler. Note the interaction 
induced diamagnetic levels currents, (b) For localized levels 
within two low energy strips. Note the horizontal logarithmic 
scale and the further suppression of PC due to interaction. 
(The visual shift between the distributions depends on the 
width of the horizontal bins). 
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the PC. This is clearly seen in Fig.^ where we consid- 
ered levels belonging to two different low energy strips 
for which the electrons were found to be localized. Such 
suppression of the PC in the localized regime by e-e in- 
teractioHS was previously seen in H-F calculations of TB 
inodelsE3, and in exact diagonalization of short range in- 
teracting electronsli^l. But exact diagonalization studies 
of TB models with Ipog range interaction show enhance- 
ment in this regimellS. We shall return to the origin of 
this discrepancy later on. 

Because the total PC in a sample is dominated by 
electrons at high levels which have the largest current 
magnitude the canonical distribution of the sample PC, 
presented in Fig.^ is clearly dominated by the single 
level distribution presented in Fi^a. The asymmetry 
in the PC direction was also analyzed in the extended 
regime in terms of the distribution of levels curvatures 
k = — ■ Here the calculations were performed by 

applying the current operator on wave functions because 
in this method one can easily calculate the levels PCs at 
flux values close to zero or to half flux quantum whick, 
for these flux values, are proportional to the curvaturestZI. 
The distributions of the curvatures were of course similar 
to the appropriate distributions in Fima. 




-1.2 I ' ' ' ' ' ' ' ' ' 1 
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FIG. 6. Comparison of tlie canonical average persistent 
current of interacting and noninteracting electrons in an en- 
semble of 150 disordered rings. The interaction effect is to 
shift the noninteracting curve towards more negative current 
values when the flux is positive. The 0o/2 periodicity is re- 
tained. However it is clear that the amplitude of the average 
current is not enhanced. 



C. Canonical average and enhancement of the single 
ring PC 

In the previous sub-sections we have presented data on 
the distribution of the canonical sample PC, single level 
PC and single level curvatures. We have seen a grad- 
ual shift in all of these distributions as a consequence 
of the e-e interactions. The sum over all electronic con- 
tributions from all samples gives the total current, and 
divided by the number of samples the average current. 
The average current as a function of flux is presented in 
Fig.^. Though the diamagnetic shift is clearly seen, the 
current amplitude is quite unchanged. This is very dif- 
ferent than the suppression found in the variouSj^gl^ 



of ID TB models in the weakly disordered regimalSt 
Nevertheless, this does not present the whole picture. For 
a particular sample a strong suppression or enhancement 
of the PC amplitude due to the e-e interactions is possi- 
ble. For example, a particular realization shown in Fig.^ 
demonstrates the highest effect of interaction on the PC 
of a single isolated ring that we found in our ensemble. 
The enhancement factor in this particular case is about 
17. Of course, even for this case, the current can not 
reach its clean ring value. 



IV. CONCLUSIONS 

Our results within the self consistent HF approximation 
for the PC carried by spinless electrons in ID disordered 
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FIG. 7. The highest enhancement factor (~ 17) of sin- 
gle ring PC due to the interaction that was found among 
our ensemble. The ring contains 42 electrons. The resulting 
PCs while calculating Veff{0) at 60 and 90 points are com- 
pared. This figure demonstrates the possibility of a huge PC 
enhancement found in the single ring experiment. 



continuous rings with long range e-e interactions can be 
summarized as follows: The PC of states in the diffusive 
and ballistic regimes are not suppressed by the e-e in- 
teractions. This is in contrast to the suppression of the 
PC found in Refs. |l|Jl| and H for the analogous TB 
models in these regimes of disorder. We found that the 
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PC amplitude for the specific samples may be altered, 
but not the amplitude of the average or typical current 
current. Further, the distribution of the PC's acquires a 
diamagnetic shift. This was confirmed by comparing the 
distributions of the individual levels currents, and also 
of the total currents, for noninteracting and interacting 
electrons. Because the total current of a sample is dom- 
inated by the extended states the average current per 
sample also shows a clear diamagnetic behavior. Our re- 
sults are significantly different from those of TB models 
where the Mott-Hubbard transition was shown to plays 
a crucial role in suppressing the PC. In the continuous 
model the Mott-Hubbard transition does not show such 
a significant role. The concept of band and filling factor 
does not seem to be relevant. In the localized regime we 
found that the PC is further suppressed by the e-e in- 
teractions. This is an artifact of the H-F approximation 
which does not describe the appearance of density corre- 
lations in this regime. An exact diagonalization of small 
ID TB rings with long-range interactioa, show that the 
PC is actually enhanced in this regimeli3. By studying 
the effect of e-e interaction on the level spacing statistics 
it is possible to gain a better understanding of the behav- 
ior of the current. The level spacing statistics show no 
sign of insulator to metal transition and this is consistent 
with the distributions found for the levels PC's. 

Introducing spin is expected to allow a stronger enhance- 
ment of the average current amplitude because of the fact 
that for spinless electrons the direct term is almost equal 
the absolute value of the exchange term, while when in- 
cluding the spin degree of freedont-|the direct term is 
about twice larger. This is crucialt3 for counteracting 
the disorder potential which leads to a higher PC am- 
plitude. This explains the important role of spin for the 
enhancement of the PC as observed_for 2D TB ringsEj 
when compared to the spinless casecj. As we have seen 
in the ID case the interactions are effective in changing 
the current in the regime of weak-disorder. One would 
expect the interaction to play an important role in the 
diffusive regime. In ID systems the states are either lo- 
calized or ballistic and therefore, the PC of the states are 
not strongly affected by the interactions. On the other 
hand, for higher dimensions where a true diffusive regime 
exists, one might expect for a stronger enhancement. 

We conclude that the effect of e-e interactions in favor 
of the enhancement of the PC in a continuous disordered 
model is stronger than in TB models. This should be 
better revealed when one considers spin and higher di- 
mensionality. We have shown that even for interacting 
spinless electrons in ID no further suppression of the PC 
due to interaction occurs. This is in contrast to the fur- 
ther suppression of the current in the weakly disordered 
TB rings that was found in Refs. |l|,E9| and U^. 
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